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Why  study  ISR? 


•  You  get  to  learn  about  many  useful  things,  in 
substantial  depth. 

-  Plasma  physics 

-  Radar 

-  Coding 

-  Electronics  (Power,  RF,  DSP) 

-  Signal  Processing 

•  But  what  if  I  probably  won’t  stay  in  this  field? 

-  See  above! 


Outline 


•  Principle  of  Pulsed  Doppler  Radar 

•  The  Doppler  spectrum  of  the  ionospheric  plasma 

•  Mathematics  of  Doppler  Processing 

•  Pulse  Compression 


The  Ubiquitous  deciBel 


The  relative  value  of  two  things,  measured  on  a  logarithmic  scale,  is  often 

expressed  in  deciBel’s  (dB) 


Example:  SNR 

Signal-to-noise  ratio  (dB)  =  10  log 
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Signal  Power 


Noise  Power 


Factor  of: 

Scientific 

Notation 
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•  dBW  =  dB  relative  to  1  watt 
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•  dBm  =  dB  relative  to  1  milliwatt 
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•  dBsm  =  dB  relative  to  1  square  meter 
of  radar  cross  section 

1,000,000 
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*  dBi  =  dB  relative  to  isotropic  radiation 

Waves  versus  Pulses 


What  do  radars  transmit? 


or  Pulses? 


by  “on-off”  action  of 
pulse  envelope 


How  many  cycles  are  in  a  typical  pulse? 

PFISR  frequency:  449  MHz] 

Typical  long-pulse  length:  480  *  5,520  cycles! 


Pulsed  Radar 

Pulse  length  100|^sec 


1  Mega-Watt 


Target 

Return 
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Inter-pulse  period 
(IPP)  i  msec 


yv 


Time 


Duty  cycle  = 


Pulse  length 


Pulse  repetition  interval 


10% 


Average  power  =  Peak  power  *  Duty  cycle 


100  kWatt 


Pulse  repetition  frequency  (PRF)  =  1/(IPP) 


1  kHz 


Continuous  wave  (CW)  radar:  Duty  cycle  =  100%  (always  on) 


*■ 
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Distance  =Time 


Range  resolution:  Set  by  pulse  length  given  in  units  of  time,  tp,  or  length,  crp 


reflected  pulse 
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Maximum  unambiguous  range:  Set  by  Inter-pulse  Period  (IPP) 


IPP  =  Interpulse  period  (s) 

PRF  =  pulse  repetition  frequence 
=  1/IPP  (Hz) 


cIPP 

u  2 


TIME- 


Velocity  =  Frequency 


Transmitted  signal:  cos(27I f  t) 


After  return  from  target:  cos 

2#„ 

2R ) 

t  + 

c  >11 

To  measure  frequency,  we  need  to  observe  signal  for  at  least  one  cycle. 

So  we  will  need  a  model  of  how  R  changes  with  time.  Assume  constant  velocity: 


Substituting: 


R-Ro  +  vot 


cos 

2k 

2v  ) 

fo  +  fo  — 

2nf0R0 

t+  ° 

\  c )  ) 

V  C  J 

~fo 


constant 


By  convention,  positive  Doppler  frequency  shift 


)  Target  and  radar  closing 
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Two  key  concepts 


Two  key  concepts: 

( - 

Distant  <^z£>Time 

R  =  cAtj  2 

Velocity  Frequency 

v  =  -fDK/2 


A  Doppler  radar  measures  backscattered  power  as  a  function  range  and  velocity. 
Velocity  is  manifested  as  a  Doppler  frequency  shift  in  the  received  signal. 


Two  key  concepts 


m/s  -60  -30  0  30  60 


A  Doppler  radar  measures  backscattered  power  as  a  function  range  and  velocity. 
Velocity  is  manifested  as  a  Doppler  frequency  shift  in  the  received  signal. 
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Concept  of  a  “Doppler  Spectrum” 


Two  key  concepts: 

( - 

Distant  <^z£>Time 

R  =  cAtj  2 

Velocity  Frequency 

v  =  -fDK/2 


Velocity  (m/s)  —  — 


If  there  is  a  distribution  of  targets  moving  at  different  velocities  (e.g.,  electrons  in  the 
ionosphere)  then  there  is  no  single  Doppler  shift  but,  rather,  a  Doppler  spectrum. 

What  is  the  Doppler  spectrum  of  the  ionosphere  at  UHF  (A  of  10  to  30  cm)? 
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Longitudinal  Modes  in  a  Thermal  Plasma 


Ion-acoustic 
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Simulated  ISR  Doppler  Spectrum 


Particle-in-cell  (PIC): 

=  —  (E(xj)  +  Vj  x  B(xj)) 
a  t  rrin 
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V-B  =  0 


Simple  rules  yield 
complex  behavior 
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ISR  Measures  a  Cut  Through  This  Surface 
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The  ISR  model 


where: 

Fe(co)=l-co  f  exp 
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From  Evans,  IEEE  Transactions,  1969 
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Incoherent  Averaging 


Normalized  ISR  spectrum  for  different  integration  times  at  1290  MHz 


0.5- 


0.5- 


0.5- 


We  are  seeking  to  estimate  the 
power  spectrum  of  a  Gaussian 
random  process.  This  requires 
that  we  sample  and  average  many 
independent  “realizations”  of  the 
process. 


Uncertainties 
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Components  of  a  Pulsed  Doppler  Radar 


s(t) 


waveform 

generator 


correlation 

receiver 


transmitter 

(amplifier) 


circulator 


Model  Fit 


p(t)  =  s(t)cos  cut 


antenna 

fill 
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low-noise 

demodulator 

amplifier  (LNA) 

p  (t)  =  a(t)  cos[  co  t  +0  (t)] 
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Plasma  density  (Ne) 

Ion  temperature  (T|) 
Electron  temperature  (Te) 
Bulk  velocity  (V|) 
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Essential  Mathematical  Operations 


Euler: 


Ae^f  =  Acos((f>)  +  jAsm(<j>) 

=  1  +  jQ 


Fourier: 


+oo 


+oo 


/(f)  = 


F{io)e^UJtdt  F(lo)  =  /  /(£)e  ^utdt 


OO 


OO 


Convolution: 


f(t)  *  .?(£)  = 


'  +  OO 


OO 


/(r)  •  g(t  -  T)dT  /(f)  *  3(t)  F(w)  •  G(w) 


Correlation: 


■  +  0O 


f(t)®g(t)=  /  /*(r)  •  g(t  +  r)dr  f(t)-kg(t)^=>F(f)**G{f) 

J  —  OO 

Wiener-Khinchine  Theorem: 


«W  ®  «*(-')  «•£/(/)•  £/*(/)  =  |  £/(/)| 
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Dirac  Delta  Function 


i  m = 


+oo9  X  =  0 


1  0, 

V 

to  ! 

S(t)  is  defined  by  the  property  that  for  all  continuous  functions 

/-boo 

5{t)f{t)dt 


—  oo 


/((- T)  =  /(()*  <5(i  -  T) 


The  Fourier  Transform  of  a  train  of  delta  functions  is  a  train  of  delta  functions. 


£  S(t  -  nT) 

n=—oo 
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Harmonic  Functions 
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Gate  function 


ect  (t/r) 


1  for  —  r/2  <  t  <  r/2 

0  otherwise 


ISR  spectrum  -< — >■  Autocorrelation  function  (ACF) 


T 


Frequency  (KHz) 


Not  surprisingly,  the  ISR  ACF  looks  like  a  sine  function... 
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Fourier  Transform  of  an  RF  Pulse 


frequency 


/o-(l/x) 


/0+(l/T) 


The  F.T.  of  a  simple  RF  pulse  is  a  sync  function  shifted  to 
the  carrier  frequency  (by  convolution  property  and 
definition  of  delta  function). 
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Fourier  Transforn 
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:h  inversely  proportional  to  the 
ind  separation  inversely 
rpulse  Period  (IPP). 
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Fourier  Transform  of  an  Long  Pulse  Train 


◄ -  T  - ► 


The  F.T.  of  an  infinite  train  of  RF  pulses  is  a  line  spectrum 
with  lines  separated  by  the  pulse  repetition  frequency. 
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Bandwidth  of  a  pulsed  signal 

Spectrum  of  receiver  output  has  sine  shape,  with  sidelobes  half  the  width  of  the 
central  lobe  and  continuously  diminishing  in  amplitude  above  and  below  main  lobe 


A 1  microsecond  pulse  has  a  null- 
to-null  bandwidth  of  the  central 
lobe  =  2  MHz 


-1  MHi  f0  + 1  MHi 


BWnn  -  2  MHz 


FREQUENCY  » 

Two  possible  bandwidth  measures: 
“null  to  null”  bandwidth  B  =  — 

nn 

T 


1 

“3dB”  bandwidth  B  ,  =  — 

^dB 

T 


Unless  otherwise  specified,  assume 
bandwidth  refers  to  3  dB  bandwidth 
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Pulse-Bandwidth  Connection 


•o 


FREQUENCY  FREQUENCY 


Shorter  pulse  ( - )  Larger  bandwidth 
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Strategies  for  radar  reception 

We  send  a  pulse  of  duration  t.  How  should  we  listen  for  the  echo? 


Convolution  of  two  rectangle  functions 


Output 


Receiver 

Input 
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Value  at  a  single  point 


J ■— I - L 


( 7,  +  7j)/2 


-to 

To  determine  range,  we  only  need  to  find  the  rising  edge  of  the 
pulse  we  sent.  So  make  T1«T2. 

But  that  means  large  receiver  bandwidth,  lots  of  noise  power,  poor 
SNR. 

Could  make  T1»T2,  then  we’re  integrating  noise  in  time  domain. 
So  how  long  should  we  close  the  switch? 
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Detection  of  a  signal  embedded  in  noise 


1  c 


1  K 
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a.  Signal  +  noise  (time  domain) 


0  100  200  300  400  500 


Sample  number 


0  0.1  0.2  0.3  0.4  0.5 

Frequency 


Exponential  pulse  buried  in  random  noise.  Sine  the  signal  and  noise  overlap 
in  both  time  and  frequency  domains,  the  best  way  to  separate  them  is  not 
obvious. 
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Most  important  consideration:  Match  the 
bandwidth  of  the  signal  you  are  looking  for 


0  10  20  30  40  50  0  0.1  0.2  0.3  0.4  0.5 

Sample  number  Frequency 


FIGURE  17-8 

Example  of  optimal  filters.  In  (a),  three  filter  kernels  are  shown,  each  of  which  is  optimal  in  some  sense.  The 
corresponding  frequency  responses  are  shown  in  (b).  The  moving  average  filter  is  designed  to  have  a 
rectangular  pulse  for  a  filter  kernel.  In  comparison,  the  filter  kernel  of  the  matched  filter  looks  like  the  signal 
being  detected.  The  Wiener  filter  is  designed  in  the  frequency  domain,  based  on  the  relative  amounts  of  signal 
and  noise  piesent  at  each  frequency. 
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The  bandwidth 


RECEIVER 

BANDWlOTH 


6.  Noise  in  receiver  output  is  proportional  to  bandwidth  of  receiver. 
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1  .  Signa|4o-noise  ratio  may  be  maximized  by  narrowing  the 

passband  of  the  IF  amplifier  to  the  point  where  only  the  bulk  of 
the  signal  energy  is  passed. 


noise  connection 


The  matched  filter  is  a  filter  whose  impulse  response, 
or  transfer  function,  is  determined  by  a  given  signal, 
in  a  way  that  will  result  in  the  maximum  attainable 
signal-to-noise  ratio  at  the  filter  output  when  both 
the  signal  and  white  noise  are  passed  through  it. 


The  optimum  bandwidth  of  the  filter,  B, 
turns  out  to  be  very  nearly  equal  to  the 
inverse  of  the  transmitted  puise  width. 


To  improve  range  resolution,  we  can 
reduce  x  (pulse  width),  but  that  means 
increasing  the  bandwidth  of  transmitted 
signal  =  More  noise... 
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Doppler  Revisited 


Transmitted  signal:  cos(27I f  t) 


After  return  from  target:  cos 

2#„ 

2R ) 

t  + 

c  >11 

To  measure  frequency,  we  need  to  observe  signal  for  at  least  one  cycle. 

So  we  will  need  a  model  of  how  R  changes  with  time.  Assume  constant  velocity: 


Substituting: 


R-Ro  +  vot 


cos 

2k 

2v  ) 

fo  +  fo  — 

2nf0R0 

t+  ° 

\  c )  ) 

V  C  J 

~fo 


constant 


By  convention,  positive  Doppler  frequency  shift 


)  Target  and  radar  closing 
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Doppler  Detection:  Intuitive  Approach 

Phasor  diagram  is  a  graphical  representation  of  a  sine  wave 
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Doppler  Detection:  Intuitive  Approach 

Closing  on  target  -  positive  Doppler  shift 

Echo  1  Echo  2  Echo  3  Echo  4  Echo  5  Echo  6  Echo  7 


Target 


Transmitted 

Received 


Echo 


Echo 


4 


5 


Echo  6  Echo  7 


Echo  1 

*  Transmitted  signal 
Echo  8 


Target’s  Doppler  frequency  shows  up 
as  a  pulse-to-pulse  shift  in  phase. 
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I  and  Q  Demodulation 


Reference  Signal 
From  Synchronizer 

T" 

* —  Detector  < — 


Reference  Signal 
Shifted  90°  in  Phase 

I 


Received 
■  Signal 
(IF) 


in-phase  (I)  channel: 

Prec(t)  COS (iOct)  —  a(t,)  COS +  U)ct)  COs(t Oct) 


=  a(t)l 


cos(</>(f )  T  2u Jet)  -(-  cos 


filter  out 


\ 


quadrature  (Q)  channel  (90°  out  of  phase): 

Prec(t)  sin(u>cf)  =  a(f)  cos (4>(t)  +  a jct)  sin(a;cf) 

/ 


=  a(t)h 


—  sin  ((f>(t)  +  2u>ct)  +  sin 


v 

filter  out 


I  and  Q  channels  together  give  the  analytic  signal 


The  fundamental  output  of  a  pulsed  Doppler  radar  is  a 
time  series  of  complex  numbers. 


1  and  Q  Demodulation  in  Free 

Transmitted  signal  Frequency  dc 

cos(2#0f) 

|uen 

jmain 

cy  C 

)omain 

-fc 

Doppler  shifted 

cos(2  n(f0  +  fD)t) 

,  c 

>  fc 

) 

-fo-fD 

COS(2  7tfj) 

cos(2  7tfDt) 

/- 

\ 

fok 

-fo-fD 

Cosine  is  even  function,  so  sign  of  fD  (and,  hence,  velocity)  is  lost. 
What  we  need  instead  is: 

expOWoO  =  cos(2  7tfDt)  +  jsm(27tfDt) 

~^D 

1 

D 

''  f0+fD 

fo 

The  analytic  signal  exp(y’2^rD?)cannot  be  measured  directly,  but  the  cos  and  sin  components 
via  mixing  with  two  oscillators  with  same  frequency  but  orthogonal  phases.  The  components 
are  called  “in  phase”  (or  /)  and  “in  quadrature”  (or  Q): 

Aexp(j2jifDt)  =  I+jQ  35 
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Example:  Doppler  Shift  of  a  Meteor  Trail 


•  Collect  N  samples  of  l(tk)  and  Q(tk)  from  a  target 

•  Compute  the  complex  FFT  of  l(tk)+jQ(tk),  and  find  the  maximum  in  the 
frequency  domain 

•  Or  compute  “phase  slope”  in  time  domain. 


Doppler  Spectra 
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Does  this  strategy  work  for  ISR? 

Typical  ion-acoustic  velocity:  3  km/s 
Doppler  shift  at  450  MHz:  1 0kHz 
Correlation  time:  I  / 1 0kHz  =  0. 1  ms 
Required  PRF  to  probe  ionosphere  (500km  range):  300  Hz 


Plasma  has  completely  decorrelated  by  the  time  we  send  the  next  pulse. 


Echo  6  Echo  7 
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Does  this  strategy  work  for  ISR? 

Typical  ion-acoustic  velocity:  3  km/s 
Doppler  shift  at  450  MHz:  1 0kHz 
Correlation  time:  I  / 1 0kHz  =  0. 1  ms 
Required  PRF  to  probe  ionosphere  (500km  range):  300  Hz 

Alternately,  the  Doppler  shift  is  well  beyond  the  max  unambiguous 
Doppler  defined  by  the  Inter-Pulse  Period  T. 
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The  ISR  target  is  “overspread” 


fd»1/T  (Doppler  changes  significantly  during  one  pulse) 


Range 


Must  sample  multiple  times  per  pulse 
Result:  Doppler  can  be  determined  from  single  pulse. 
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ISR  Receiver:  Doppler  filter  bank 


Doppler  Filter  Bank  Doppler  Range 


Choose  fL  and  fH  such  that 
n_<kr'fdmax  and  fH'>fo"*"Tdmax 


Practical  Problem:  It  is  hard  to  make  narrow  band  (High  Q)  RF  filters: 
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ISR  Receiver:  I  and  Q  plus  decorrelation 


We  have  time  series  of  V(t)  =l(t)  +  jQ(t),  how  do  I  compute  the  Doppler  spectrum? 


Estimate  the  autocorrelation 
function  (ACF)  by  computing  products 
of  complex  voltages 


p(T)  = 


(“lag  products”) 

A _ r 

(V(t)V‘(t  +  T)) 

Power  spectrum  is  Fourier 
Transform  of  the  ACF 
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Pulse  compression  and  matched 

"If  you  know  what  you're  looking  for,  it's  easier  to  find." 


100  pixels 
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Range  detection:  revisited 


T  =  Pulse  length 
B  =  Bandwidth 


For  high  range  resolution  we  want  short  pulse  O  large 
bandwidth 


For  high  SNR  we  want  long  pulse  O  small  bandwidth 


Long  pulse  also  uses  a  lot  of  the  duty  cycle,  can’t  listen  as 
long,  affects  maximum  range 

The  Goal  of  pulse  compression  is  to  increase  the 
bandwidth  (equivalent  to  increasing  the  range  resolution) 
while  retaining  large  pulse  energy. 


Linear  Frequency  Modulation  (LFM  or 
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Matched  filter  detection  of  a  chirp 


RECEIVED 


FREQUENCY 

- ► 


Since  trailing  portions  of  echo  take  less  time  to  pass  through 
filter,  successive  portions  tend  to  bunch  up:  Amplitude  of 
pulse  is  increased  and  width  is  decreased. 


b  +  A 


Ilf  I 


B  |8|7|6|5|d|3|2|l~| 

A  |8|7|6|5|4|3|2|l' 


Echoes  from  closely  spaced  targets,  A  and  B,  are  merged  but, 
because  of  coding,  separate  in  output  of  filter. 
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Barker  codes 


+  +  + 


correlator  output 


+ 


+  + 


+  + 


+  +  + 


+ 


-1+1=0 
11+1=1 
1  +  1-1-1=0 
1+1+1 +1+1=5 


TABLE  6.2  All  Known 
Binary  Barker  Codes 


Code  Length 

Code 

2 

11  or  10 

3 

110 

4 

1110  or  1101 

5 

11101 

7 

1110010 

11 

11100010010 

13 

1111100110101 

‘-'f  p ^ -i'  ^ ^  ^ ^ 
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Matched  filtering  of  Barker  Code 


TIMEO 


TIME  1 


TIME  2 


TIME  3 


TIME  4 


TIME  5 


TIME  6 


DELAY  LINE 


+  +  - 
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0UTPU1 


-►  0 
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-►  +3 


3  bit  code 


TIME 


TABLE  6.2  All  Known 
Binary  Barker  Codes 


Code  Length 

Code 

2 

11  or  10 

3 

110 

4 

1110  or  1101 

5 

11101 

7 

1110010 

11 

11100010010 

13 

1111100110101 

1 5.  Step-by-step  progress  of  o  3-digit  binary  phase  modulated  pulse  through  a  tapped  delay  line. 
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